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量子 鲸鱼 优化 算法 求解 作业 车 间 调度 问题 
间 ” 旭 ， 叶 春明 ， 姚 远 运 


(上 海 理工 大 学 管理 学 院 ， 上 海 200093) 


摘 要 : 为 了 克服 基本 鲸鱼 优化 算法 (WOA) 在 解决 作业 车 间 调 度 问 题 时 存在 收敛 精度 低 、 容 易 陷 入 局 部 最 优 的 缺陷 ， 
利用 量子 计算 与 优化 思想 提出 了 一 种 量子 综 鱼 优化 算法 (QWOA)， 并 对 其 进行 了 计算 复杂 度 分 析 、 全 局 收敛 性 证 明 及 
仿真 实验 。 通 过 对 11 个 作业 车 间 调度 问题 基准 算 例 的 仿真 实验 发 现 ， 与 基本 镶 鱼 优化 算法 (WOA)、 布 谷 鸟 搜索 算法 
(CS)、 灰 狼 优 化 算法 (GWO) 相 比 ，QWOA 算法 在 最 小 值 、 平 均值 、 寻 优 成 功率 等 方面 具有 较 优 结果 。 研 究 表明 ， 
量子 鲸鱼 优化 算法 在 解决 作业 车 间 调度 问题 时 ， 具 有 更 高 的 收 系 精度 和 更 好 的 全 局 搜索 能 力 ， 且 能 够 跳出 局 部 最 优 。 
关键 词 : 鲸鱼 优化 算法 ; 量子 计算 与 优化 ; 作业 车 间 调度 ; 收 你 性 证 明 ; 混合 算法 
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Solving Job-Shop scheduling problem by quantum whale optimization algorithm 


Yan Xu, Ye Chunming, Yao Yuanyuan 
(Business School, University of Shanghai for Science and Technology, Shanghai 200093, China) 


Abstract: To overcome whale optimization algorithm’s disadvantages of poor convergence and being easily trapped in local 
optima in solving job-shop scheduling problem, this paper proposed a quantum whale optimization algorithm (QWOA) based 


on the idea of quantum computation. Then it provided the computational complexity analysis, global convergence proof and 


simulation experiments of QWOA. Based on a set of 11 JSP benchmark instances, the simulation experiments show that QWOA 
© has better results than whale optimization algorithm (WOA) , cuckoo search (CS) and grey wolf optimizer (GWO) , in the 
Se minimum value, average value and success rate. Finally, this paper concludes that QWOA has the merits of higher convergence 


accuracy, higher local optima avoidance and better exploration. 


Key Words: whale optimization algorithm; quantum computing and optimization; job-shop scheduling problem; convergence 


proof; hybrid algorithm 
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失 容 易 陷入 局 部 最 优 、 收 敛 速度 慢 、 收 敛 精度 低 等 方面 的 不 足 ， 
有 具体 涉及 种 群 初始 化 方法 的 选择 与 优化 、 种 群 多 样 性 的 保持 与 

作业 车 间 调 度 问题 (Job-Shop scheduling problem，JSP) 是 提升 、 局 部 开发 能 力 及 全 局 搜索 能 力 的 平衡 与 提升 、 可 行 解 的 
最 经 典 的 离散 制造 系统 调度 问题 。 作 为 一 种 最 困难 的 组 合 优化 。” ”空间 扰动 与 遍历 等 方面 ; 如 , Jangir 等 人 四 在 原 有 算法 的 基础 上 
问题 ，JSP 已 经 被 证 明 是 一 种 NP-hard 问题 ， 在 多 项 式 时 间 提出 了 自 适 应 鲸鱼 优化 算法 (AWOA)，Ling 等 人 外 提出 了 
找 不 到 一 个 算法 精确 地 求解 其 最 优 解 ， 该 问题 的 求解 算法 和 LWOA 算法 , 钟 明 辉 等 人 09 提 出 了 一 种 随机 调整 控制 参数 的 高 


0 引言 


湾 王 


已 成 为 车 间 调 度 领 域 的 热点 由。 因此 ， 本 文 将 采用 新 兴 的 鲸鱼 效 的 鲸鱼 优化 算法 (EWOA), 刘 竹 松 等 人 0 提出 了 正 余 混 沌 双 
优化 算法 解决 作业 车 间 调 度 问 题 。 弦 饼 鱼 优化 算法 (CSCWOA)。 综 上 所 述 , 馈 鱼 优化 算法 多 用 于 
鲸鱼 优化 算法 (whale optimization algorithm，WOA)， 连续 优化 问题 求解 ， 且 尚未 被 应 用 于 解决 作业 车 间 调 度 问题 。 


Mirjalili 等 人 所 于 2016 年 提出 ， 其 仿生 学 原理 源 自 于 座 头 鲸 的 因此， 本 文 将 借鉴 量子 计算 与 优化 的 思想 对 鲸鱼 优化 算法 进行 
狩猎 行为 。 目 前， 对 于 鲸鱼 优化 算法 的 研究 尚 处 于 起 步 阶段 ， 改进 ， 并 将 其 用 于 解决 离散 的 作业 车 间 调 度 问题 。 


其 正在 被 逐渐 应 用 到 更 多 的 领域 ， 如 资源 调度 领域 Bl、 选 址 与 
1 业 车 间 调 度 问 题 的 数学 描述 
路 径 规划 领域 9、 工 业 设计 领域 [7 等 。 在 算法 的 改进 方面 , 其 作业 车 间 调度 问题 的 数学 描述 
主要 研究 目标 为 对 算法 的 收敛 精度 和 收敛 速度 的 提升 ， 克 服 算 作业 车 间 调度 问题 (JSP) 一 般 被 描述 为 : 在 给 定 每 个 工件 
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的 加 工 工艺 、 每 个 工件 使 用 机 器 的 序列 及 每 个 工件 每 道 工 序 的 
加 工时 间 的 情况 下 ， 安 排 一 种 工件 的 加 工 顺序 ， 使 得 个 待 加 
工 的 工件 在 m 台 机 器 上 进行 加 工 ， 满 足 某 种 性 能 指标 最 优化 。 
本 文 所 使 用 的 性 能 指标 为 “最 小 化 最 大 完工 时 间 ” 且 该 问题 需 
要 满足 如 下 假设 及 约束 : 
a) 机 器 不 发 生 故 障 ， 从 零 时刻 开 始 加 工 工件 。 
b) 同一 台 机 器 在 同一 时 刻 只 能 加 工 一 个 工件 ; 同一 工件 同 
时 刻 只 能 在 一 台 机 器 上 加 工 。 
c) 工 件 的 加 工 工 序 严格 遵守 加 工时 间 , 一 旦 开始 不 得 中 断 ; 
工序 间 不 允许 抢占 加 工 。 
d) 每 个 工件 的 工艺 路 线 固 定 , 即 工件 使 用 机 器 的 序列 给 定 
每 个 工件 在 工艺 路 线 中 的 每 个 工序 的 加 工时 间 固 定 。 但 ， 不 同 
工件 的 工序 间 无 顺序 约束 。 
eg nn 术 有 线性 规划 模型 
(LP) 、 整 数 规划 模型 (IP〉 描述 、 析 取 图 模型 描述 。 现 
ee 


min max{cx} 


A 


一 


Es” 有 

St. Cx -tr +M(l—anm)>c,, Q2) 

RC +M(1- i) ti G3) 

cx >0 (4) 

ling Xi = 0,1 (5) 

i, j=1,2,..n ; hk =1,2,...,m (0) 

其 中 : 如 和 cy 表示 工件 ;在 机 器 上 的 加 工时 间 和 完工 时 间 。 
式 (5) 中 ，aww =1 表 示 机 器 有 先 于 机 器 加 工 工件 i; xz =1 
表示 工件 ; 先 于 工件 j 在 机 器 上 加 工 ; 反之 , 二 者 省 为 0。jy 
一 个 大 数 ， 则 式 〈2 ) 表示 工序 的 前 后 约束 关系 ; 式 (3) 表 


示 工 序 的 非 阻 塞 约 束 关系 ; 
最 大 完工 时 间 的 最 小 化 0 。 
2 ”量子 鲸鱼 优化 算法 
2.1 基本 鲸鱼 优化 算法 (WOA) 

WOA 算法 的 仿生 学 原理 受 启发 于 座 头 鲸 独 特 的 砚 食 行为 
(bubble-net feeding method)。 其 原理 为 ， 座 头 鲸 发 现 融 食 目标 
后 ， 沿 形似 螺旋 的 路 笃 ， 通 过 制造 气泡 网 的 方式 进行 员 食 ， 上 县 
体 行为 包含 “气泡 网 攻击 (bubble-net attacking method )”“ 收 
缩 包围 (encircling prey)”“ 随 机 搜寻 (Search for prey)” 等 三 
种 主要 方式 。 其 中 气泡 网 攻击 捕食 猎物 的 方式 又 可 依据 概率 分 
为 “螺旋 式 位 置 更 新 (spiral updating position )” 和 “缩小 包围 
则 座 头 鲸 的 


式 (1) 为 目标 函数 ， 表 示 性 能 指标 


距离 (shrinking encircling mechanism )” 两 个 步骤 。 
员 食 行为 可 进行 如 下 数学 描述 
2.1.1 气泡 网 攻击 

鲸鱼 以 群体 方式 砚 食 ， 砚 食 过 程 中 ， 鲸 鱼 个 体 根 据 随 机 概 
率 值 的 大 小 ， 对 “螺旋 式 位 置 更 新 (spiral updating position: )” 
和 “缩小 包围 距离 〈shrinking encircling mechanism )” 这 两 种 行 


为 进行 选择 。 其 中 ,“ 缩 小 包围 距离 (shrinking encircling 
mechanism)” 与 2.1.2 节 表 述 相 同 ,“ 螺 旋 式 位 置 更 新 (spiral 
updating position )” 的 数学 模型 可 描述 为 


X(1+1)=D’.e" .cos(2x1)+ X(t) 0 


D=|x (0)-X() (8) 


其 中 :“. ”表示 点 乘 ,“ 


”表示 取 绝 对 值 ，1 表示 当前 的 迭代 


次 数 ， 向 量 匀 (1) 表示 座 头 馈 的 当前 位 


， 向 量 X(1+1) 表示 座 


融 食 位 


头 鲸 办 代 更 新 后 的 位 置 ， 向 量 这;(r) 表示 鲸 群 中 的 最 佳 


置 ， 参 数 b 为 常数 用 来 定义 对 数 螺旋 的 形状 ， 参 数 ] 为 [-1,1] 之 
间 的 随机 数 ， 两 者 共同 控制 鲸鱼 个 体 的 螺旋 式 位 置 更 新 方式 。 

当 1= -1 时 鲸鱼 个 体 距 离 邵 食 目 标 最 近 ， 当 7 =1 时 鲸鱼 个 体 距 
见 员 


标 最 远 。 
在 气泡 网 攻击 过 程 中 ， 将 鲸鱼 个 体 选 择 位 置 更 新 方式 的 行 
为 概率 设 为 p 。 则 气泡 网 攻击 (Bubble-net attacking method ) 过 
程 可 以 完整 描述 为 


ed 


X(t)-A.D, 
D’'.e” ‘cos(2x1)+ xX’ (1) 


p<p” 
p>p" 


X(t1+1)= (9) 


2.1.2 收缩 包围 (encircling prey) 
根据 2.1.1 节 描述 ， 当 P< 天 时 〈 产 取 05)， 鲸 鱼 个 体 更 新 


自身 在 空间 中 的 位 置 ， 向 较 优 的 位 置 移动 ， 以 实现 对 猎物 的 收 
缩 包 围 。 其 数学 模型 描述 如 下 : 
5=|c-x()-X() (10) 
X(t+1)=X (1)-A.D (11) 
A=2a:n-a (12) 
C=2n (13) 
21 
4S (14) 


其 中 : 区, 为 最 大 连 代 次 数 ，a 的 值 随 着 迭代 次 数 的 增加 由 2 
向 0 线性 递减 ，# 和 万 为 (0,1) 的 随机 向 量 ， 系数 向 量 忆 入 
控制 鲸鱼 个 体 的 游 走 方式 。 需 要 指出 的 是 ， 鲸 鱼 个 体 采 取 此 种 


位 置 更 新 的 前 提 是 4 


小 于 工时 母 鱼 个 体 又 (1) 


小 于 1。 即 , 当 |4 


向 统 群 中 的 当前 最 优 位 置 (1) 靠近 。 


2.1.3 随机 搜寻 (search for prey) 


在 2.1.2 节 中 ， 


式 (11) 进行 位 置 更 新 的 前 提 是 加 | 小 于 1。 


不 小 于 1 时 , 采取 如 下 随机 方式 进行 位 置 更 新 , 并 称 之 


则 当 |4 
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为 随机 搜寻 (search for prey)。 数 学 模型 描述 如 下 : 


X(t+1)= Xn - A: Da (15) 


Dr =|C: Xm -X() (16) 


其 中 ， 元 ,表示 鲸鱼 群 中 随机 的 其 他 鲸鱼 个 体 。 式 (15) (16) 


表示 当 [4 不 小 于 1 时， 鲸鱼 个 体 随机 向 其 他 个 体 游 走 靠 近 


2.2 ”量子 鲸鱼 优化 算法 (QWOA) 
2.2.1 量子 计算 与 量子 优化 
量子 计算 (quantum computation，QC ) 


上 Feynman 首先 于 
1982 年 提出 , 现 已 成 为 国内 外 学 者 关注 的 前 沿 学 科 。 在 量子 计 
算 中 ,信息 的 基本 存储 单位 称 为 量子 比特 (quantum bit)， 有 别 
于 经 典 比 特 ， 量 子 比 特 可 以 连续 、 随 机 地 存储 于 两 个 极 化 状态 
的 任意 亚 加 态 ， 按 类 型 可 分 为 单 量子 比特 、 双 量子 比特 、 多 量 
子 比 特 。 其 中 单 量子 比特 的 登 加 状态 可 以 表示 为 
DO)=awlO+O| 1) (17) 


及 


Ilo 3 本 


中 : |0) 及 上 1 为 量子 计算 中 微观 粒子 的 两 种 基本 状态 的 狄 拉 
克 记 号 表示 ，w 及 1 表示 量子 态 的 概率 幅 ， 为 一 对 复数 。 
量子 计算 中 对 量子 状态 进行 逻辑 变换 的 量子 装置 称 为 量子 
逻辑 门 ， 是 实现 量子 计算 的 基础 。 常 用 的 量子 门 有 “Hadamard 
门 ”“Pauli-X 门 ”“Pauli-Y | 门 ”“Pauli-Z 门 ”“ 相 位 门 ”“m/8 门 ” 
“量子 旋转 门 ” 等 。 其 中 量子 旋转 门 可 表示 为 

“a 
子 状态 的 更 新 过 程 可 表示 为 


Ho 


其 中 。 6 为 有 转角 度 ，| | 为 更 新 后 的 概率 1 


开 


上 
其 量子 


2.2.2 量子 鲸鱼 优化 算法 

于 基本 鲸鱼 优化 算法 存在 收敛 精度 低 、 容 易 陷 入 局 部 最 
缺点 ,因此 采用 量子 计算 与 量子 优化 的 思想 对 其 进行 改进 ， 
更 新 方 
己 知 


优 和 
将 量子 旋转 门 操作 加 入 基本 鲸鱼 优化 算法 中 ， 对 其 位 
式 进行 改进 ， 以 提升 基本 算法 的 种 群 多 样 性 和 收敛 精度 。 


丰 


叉 (1) 为 馈 鱼 在 搜索 空间 中 的 位 


, 设 OX (1) 为 鲸鱼 经 过 量子 旋 


转 门 操作 后 的 位 置 。 
式 (18) 可 灸 


可 以 简化 为 


知 ，5=VL-o2， 则 量子 旋转 门 的 更 新 过 程 


a =|acos(0)- Mi-e ‘sin(9) (21) 


相对 式 (20) 而 言 ， 式 (21) 的 量子 旋转 门 操作 方式 被 称 
为 单 链 编码 量子 旋转 门 操作 5。 鲸鱼 个 体 经 过 单 链 编码 量子 旋 


we 02) 


当 鲸 鱼 个 体 经 过 量子 旋转 门 更 新 后 ， 将 其 所 对 应 角 
与 旋转 之 前 解 的 质量 进行 对 比 ， 并 选择 其 中 较 优 者 。 灿 
鱼 优化 算法 简称 之 为 QWOA(quantum whales optimization 
algorithm)， 其 具体 步骤 可 分 为 : 

0) 设置 迭代 次 数 、 种 群 规模 、 问 题 参 数 。 

b) 利 用 随机 初始 化 方法 生成 初始 种 群 。 

9 计算 鲸鱼 个 体 的 适应 度 值 ， 找 出 种 群 中 的 当前 全 局 最 优 
解 的 位 置 


| 
= 
| 妆 


Ey 

二 
内 
| 
时 


dd) 根据 p 及 4 的 取 值 情况 对 鲸鱼 个 体 的 每 个 维度 的 变量 


进行 位 置 更 新 。 
日 对 鲸鱼 个 体 各 维度 的 变量 进行 量子 旋转 门 操作 。 
们 计算 d) 与 9) 的 更 新 效果 , 并 对 二 者 进行 选择 , 择优 保 
多 在 最 大 和 迭代 次 数 内 ， 对 鲸鱼 群体 重复 c)- 人 操作 。 
h) 输出 求解 结果 。 
至 此 ， 量 子 鲸鱼 优化 算法 的 流程 可 完整 描述 如 图 1 所 示 。 


可 


Stepl: 
设置 算法 种 群 规模 、 
迭代 次 数 及 问题 参数 


Step2: 
进行 种 群 初始 化 操作 


是 否 小 于 最 有 Step8: 
大 迁 代 次 数 一 输出 最 终 
运算 结果 
是 
YY 
Step3: 
找 出 当前 全 局 最 优 解 
Step4 
Step7 p>=p*? 是 
否 
按 式 (15) 进行 按 式 (11) 进行 | | 按 式 〈7) 进行 
位 置 更 新 位 置 更 新 螺旋 式 位 置 更 新 
Step5: 
按 式 〈22) 进行 量子 旋转 门 操作 更 新 位 置 
Step6: 


法 流程 


图 1 量子 鲸鱼 优化 


3 ”计算 复杂 度 分 析 及 收敛 性 证 明 


3.1 计算 复杂 度 分 析 
利用 频 度 估计 的 方法 对 量子 鲸鱼 优化 算法 进行 计算 时 间 复 
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杂 度 和 空间 复杂 分 析 。 设 算法 的 最 大 和 迭代 次 数 为 7 ， 种 群 规模 
为 POP，, 变量 的 维度 为 pjy 。 由 2.2.2 节 可 知 , Step4- Step6 为 
算法 共同 的 原 操 作 ， 为 处 于 最 内 层 循环 的 循环 体 ， 语 句 的 执行 
频 度 均 为 f(n)=T.POP.DIM ， 则 算法 的 渐进 时 间 复 杂 度 
T(n) 
受到 种 群 规 模 POP 和 变量 维度 DIM 的 影响 ， 则 其 空间 复杂 度 


=0O(f(n))=0(7T.POP.DIM) 在 算法 中 , 可 变 的 储存 空间 


可 表示 为 s(n)=0O(f(n))=0O(POP. DIM) 


3.2 ”算法 收敛 性 证 明 05 
利用 概率 测度 法 对 QWOA 算法 进行 全 局 收敛 性 证 明 。 根 


据 全 局 收敛 性 准则 及 定理 ， 要 证 明 算 法 为 全 局 收敛 算法 ， 则 算 
法 需要 满足 如 下 两 个 假设 : 
假设 1 f(D(z,é))<f(z)， 并 且 如 果 ges， 则 

5S f(D(60)s7(0): 
< 。 其 中 : / 为 最 小 化 问题 的 目标 函数 ，D 为 可 以 保证 产生 较 优 最 
忆 新 解 的 算 子 或 函数 。z 为 解 空间 Rr 的 子 集 $ 的 一 个 点 , 能 够 使 
过。 得 函数 的 值 最 小 化 或 在 $ 上 产生 可 接受 的 函数 值 的 下 确 界 。 < 
LC ”为 算法 在 样本 空间 上 生成 的 向 量 。 此 假设 的 意义 在 于 保证 优化 
中 ”算法 的 正确 运行 ， 以 使 算法 的 解 序列 收敛 于 函数 f 在 s 上 的 下 
ee 
于 确 界 , 。 
© 假设 2 ”对 于 5 的 任意 Borel 子 集 4 ， 若 其 概率 测度 
| 
> Y[ 久 >0,， 则 有 TI(4-A[4)=0 
泛 。 其 中 : 44 [4] 为 由 分 布 产生 的 对 4 的 概率 测度 。 此 假设 的 意 
它 。 义 在 于 ， 在 测度 大 于 0 的 情况 下 ， 算 法 经 无 穷 多 次 迭代 后 ， 不 
-三 。 会 错过 解 空间 8 的 任意 Borel 子 集 4， 
© 按 上 述 两 假设 要 求 对 QWOA 算法 的 全 局 收敛 性 进行 证 明 。 


引 理 1 QWOA 满足 假设 1。 
证 明 1 根据 2.2.2 节 描 述 ， 可 将 函数 刀 在 QWOA 算法 中 
的 描述 定义 为 


G, f(s(X,))=/f(G6,) 
gs(Xi,), f(s(Xi))<f(G,) 

D(G.,X,,)= Ga C3) 
g(X)s 下 可 和 JJj<7(G) 


纹 鱼 群 当前 全 局 最 优 位 置 解 的 计算 与 选择 ， 


其 中 : 函数 p 表示 对 


连续 应 用 


函数 g 对 应 量子 旋转 操作 函数 ，g (X,, ) 表 示 鲸 鱼 个 体 


i 的 第 1 次 更 新 后 的 位 置 ，G, 为 当前 的 全 局 最 优 解 的 位 置 。 按 

照 本 文 算法 的 定义 ， 可 知 G 所 对 应 的 适应 度 值 是 单调 不 增 的 ， 
逐渐 向 解 空间 的 下 确 界 收敛 ， 因 此 引 理 1 得 证 。 
引 理 2 QWOA 满足 假设 2。 


g(xX,,,) = [Xi :cos(O) - .1- X,j, ‘sin(0)| (24) 


其 中 : X 为 任意 迭代 步 1 ,第 i 只 鲸鱼 个 体 的 第 ) 维度 变量 的 


位 置 ，g(X ，) 、Xij, 的 取 值 范围 为 [0,1]。 则 由 式 (24) 可 得 : 
0=arccos(+g (Xj,))-o (25) 
其 中 : co =arccos(X, ,,) 的 取 值 范围 为 [0 


7] 欲 使 式 (25) 成 


立 ， 则 g 的 取 值 范 填 


dT 
ea 
IE 


辽 ]， 设 X 为 Rw 第 7 维 所 ， 


g(X,,,) 为 Rw 第 j 维 不 同 于 X,,, 的 任意 点 。 可 知 当量 子 旋转 


角度 落 在 0 范围 内 时 ， Xi 到达 sg( 


XX ) 的 概率 大 于 0。 


定义 44, 为 对 应 于 鲸鱼 个 体 N 维 变量 的 均匀 分 布 ， 由 上 述 
可 知 ， 对 8 的 任意 Borel 子 集 4， 当 Y[4]>0 ， 可 得 


0<w,[4] <1 (26) 
进而 可 得 ， 由 六 产生 的 对 4 的 概率 测度 : 
4[4]=1 -I1( -4 [4]) (27) 
通过 式 (26) (27) 得 到 
I[(-4[A)=0 Cg) 


至 此 ， 引 理 2 得 证 。 
所 以 ， 根 据 证 明 1 及 证 明 2 可 证 QWOA 算法 满足 假设 1 
及 假设 2， 为 全 局 收敛 算法 。 


4 ”仿真 实验 与 分 析 


为 了 确定 和 检验 鲸鱼 优化 算法 解决 作业 车 间 调 度 问 题 的 性 
能 , 本 文选 取 作 业 车 间 调 度 标准 测试 问题 库 中 的 11 个 测试 问题 
进行 求解 并 与 文献 9 中 的 其 他 算法 作出 比较 。 算 法 的 运算 环境 
为 Windows10 专业 版 32 位 操作 系统 ， Intel(R) Core(TMD i7- 
2600 CPU @ 3.40GHz 处 理 器 ，4GRAM。 算 法 的 编程 实现 软件 
为 Matlab2015b 。 
4.1 WOA 与 灰 狼 优化 算法 和 布谷 鸟 搜 索 算 法 的 比较 

灰 狼 优化 算法 07 (grey wolfoptimizer，GWO) 由 Mirjalili 
等 人 于 2014 年 提出 ; 布谷 鸟 搜索 算法 08 (cuckoo search，CS ) 
Yang 等 于 2009 年 提出 ; 二 者 皆 为 新 兴 的 智能 算法 。 为 了 更 
准确 地 与 GWO 和 CS 算法 就 JSP 问题 进行 比较 ， 将 WOA 算 
去 在 求解 JSP 问题 时 的 参数 设置 与 文献 9 保持 一 致 ， 即 种 群 规 
模 为 30， 最 大 迭代 次 数 为 500， 算 法 独立 运行 次 数 为 30。 得 到 
如 表 1 所 示 结 果 。 


表 1 的 结果 可 知 ， 在 种 群 规模 、 和 迭代 次 数 及 算法 运算 次 
数 相同 的 条 件 下 ，WOA 算法 与 GWO 算 和 CS 算法 在 FT06、 


二 < 
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LA01、LA06、LA11 等 问题 中 都 能 收敛 到 当前 已 知 最 优 解 ( 部 
分 甘 特 图 如 图 2 所 示 )， 但 在 寻 优 成 功率 上 WOA 算法 要 劣 于 
GWO 算法 和 CS 算法 。 


在 最 小 值 方面 ，WOA 算法 优 于 CS 算法 ， 但 劣 于 GWO 算 
法 。 在 平均 值 方面 ，WOA 算法 劣 于 GWO 算法 ， 少 部 分 优 于 
CS 算法 。 


表 1 基本 鲸鱼 优化 算法 与 灰 狼 优化 算法 和 布谷 鸟 搜索 算法 的 比较 


实例 名 称 WOA GWO CS 
己 知 最 优 

(问题 规模 ) 最 小 值 ”平均 值 ” 寻 优 成 功率 ”最 小 值 ” 平 均值” 和 寻 优 成 功率 ”最 小 值 ” 平 均值” 寻 优 成 功率 
FT06 (6x6) 55 55 57.8 16.7% 55 57.6 20% SS 55.0 96.7% 
FT10 (10x10) 930 1065 1155.2 0% 989 1061.8 0% 1075 1131.6 0% 
FT20 (20x5) 1165 1334 1467 0% 1274 -1378.4 0% 1368 1429.2 0% 
LAO1 (10x5) 666 666 689.8 30% 666 673.7 56.7% 666 667.0 86.7% 
LA06 (15x5) 926 926 933.8 53.3% 926 926.6 96.7% 926 926 100% 
LA11 (20x5) 1222 1222 1241.0 40% 1222 1222.1 93.3% 1222 1226.8 56.7% 
LA16 (10x10) 945 994 1080.9 0% 979 1026.4 0% 1015 -1043.8 0% 
LA21 (15x10) 1046 1256 1337.6 0% 1167 1235.3 0% 1267 1331.5 0% 
LA26 (20x10) 1218 1510 1598.7 0% 1345 1459.3 0% 1487 1614.6 0% 
LA31 (30x10) 1784 1978 2161.1 0% 1871 1979.8 0% 2102 2198.8 0% 
LA36 (15x15) 1268 1555 1658.5 0% 1393 1501.7 0% 1600 1655 0% 


0 10 20 30 40 50 55 


时 间 (min) 


标准 测试 问题 LA01- 甘 特 


0 100 200 300 400 500 600 666 
时 间 (min) 


图 2 部 分 标准 


E 测 试问 题 的 甘 特 图 


上 述 结果 说 明 ，WOA 算法 虽然 能 够 成 功 运 用 于 在 作业 车 
间 调 度 问题 中 ， 但 与 GWO 与 CS 算法 相 比 较 并 无 明显 优势 ， 
仍然 存在 较 大 改进 空间 。 因此, 本 文采 用 改进 后 的 QWOA 算法 
继续 对 上 述 问 题 进行 求解 。 
4.2 QWOA 与 WOA 的 比较 

利用 量子 鲸鱼 优化 算法 对 上 述 问题 进行 求解 ， 将 算法 的 种 
群 规模 设置 为 30， 最 大 返 代 次 数 设置 为 500， 量 子 旋转 的 角度 
设置 为 rand*2* x 。 运 算 结果 的 对 比如 表 2 所 示 。 

表 2 中 ， 最 小 值 提升 、 平 均值 提升 、 寻 优 成 功率 提升 的 计 
算 方法 分 别 为 : (WOA 最 小 值 -QWOA 最 小 值 )/ WOA 最 小 值 、 

(WOA 平均 值 -QWOA 平均 值 ) /WOA 平均 值 、QWOA 寻 优 

成 功率 -WOA 寻 优 成 功率 。 

通过 QWOA 与 WOA 算法 就 求解 作业 车 间 调 度 问 题 的 比 
较 可 知 ， 量 子 鲸鱼 优化 算法 对 算法 的 寻 优 成 功率 方面 有 很 大 的 
提升 , 在 FT06、LA01、LA06、LA11 等 问题 上 的 寻 优 成 功率 的 


值 均 在 83% 以 上 , 其 中 LA01、LA06、LA11 的 寻 优 成 功率 分 别 
高 达 93.3%、100%、100%, 寻 优 成 功率 的 平均 提升 率 为 59.2%; 
且 在 LA31 问题 中 ,收敛 到 了 WOA、GWO、CS 算法 所 无 法 收 
敛 到 的 已 知 最 优 解 ， 同 时 在 LA16 问题 中 也 近乎 收敛 到 已 知 最 
优 解 。 

而 且 ，QWOA 对 最 小 值 及 平均 值 的 提升 效果 也 十 分 显著 ， 
取 值 明显 优 于 WOA、GWO、CS 算法 , 对 WOA 算法 的 最 小 值 
及 平均 值 的 平均 提升 率 分 别 为 10.5%、9.7%。 

同时 ， 图 3 中 QWOA 与 WOA、CS、GWO 等 算法 寻 优 曲 
线 的 特征 对 比 ， 说 明 量子 鲸鱼 优化 算法 在 求解 作业 车 间 调 度 问 
题 的 最 优 解 时 , 收敛 到 最 优 值 所 需要 的 迭代 次 数 更 少 ;收敛 的 精 
度 更 高 ， 跳 出 局 部 最 优 进行 全 局 搜索 的 能 力 更 强 。 

结合 算法 寻 优 的 方差 情况 ( 受 篇 幅 限 制 未 予 展 示 ), 总 体 而 
言 ， 利 用 量子 计算 与 量子 优化 的 思想 对 基本 鲸鱼 优化 算法 进行 
改进 后 ， 量 子 鲸鱼 优化 算法 在 求解 作业 车 间 调 度 优 化 问题 时 有 具 
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人 
_ _ChinaXiv 合 作 期 十 
间 旭 ， 等 : 量子 鲸鱼 优化 算法 求解 作业 车 间 调度 问题 
有 更 高 的 收敛 精度 、 更 少 的 收敛 次 数 和 更 强 的 全 局 搜索 能 力 和 ” ” 和 鲁 棒 性 。 
表 2 量子 鲸鱼 优化 算法 与 基本 鲸鱼 优化 算法 的 比较 
QWOA WOA 改进 提升 效果 
实例 名 称 i 
oe 最 优 寻 优 寻 优 寻 
(问题 规模 ) 最 小 值 ”平均 值 最 小 值 ”平均 值 最 小 值 ”平均 值 
FT06 (6x6) SS SS SS.3 83.3% SS 57.8 16.7% 0% 4.3% 66.6% 
FT10 (10x10) 930 966 1007.2 0% 1065 1155.2 0% 9.3% 12.8% 0.0% 
FT20 (20x5) 1165 1207 1252:] 0% 1334 1467 0% 9.5% 14.6% 0.0% 
LAO1 (10x5) 666 666 667.5 93.3% 666 689.8 30% 0% 3.2% 63.3% 
LA06 (15x5) 926 926 926 100% 926 933.8 53.3% 0% 0.8% 46.7% 
LA11 (20x5) 1222 1222 1222. 100% 1222 1241.0 40% 0% 1.5% 60.0% 
LA16 (10x10) 945 946 994.3 0% 994 1080.9 0% 4.8% 8.0% 0.0% 
LA21 (15x10) 1046 1110 1146.2 0% 1256 1337.7 0% 11.6% 14.3% 0.0% 
LA26 (20x10) 1218 1269 1331.1 0% 1510 1598.7 0% 16% 16.7% 0.0% 
LA31 (30x10) 1784 1784 1846.2 6.7% 1978 2161.1 0% 9.8% 14.6% 6.7% 
LA36 (15x15) 1268 1357 1403.8 0% 1555 1658.5 0% 12.7% 15.4% 0.0% 
标准 测试 问题 FT06- 寻 优 曲 线 标准 测试 问题 FT20- 寻 优 曲线 
75 厂 一 1900 
| E 
| 二 1800 
7ol 
! 1700 
i i 
= > 1600 
要 袁 1500 
对 "| 
1400 
60 
1300 
55 1200 
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5 


鲸鱼 


续 优 化 问题 ,在 作业 车 间 调 
和 。 本 文 验 训 


较 少 应 | 


结束 语 


优化 算法 (WOA) 是 新 兴 的 群 智能 算法 ， 多 应 用 于 
度 (JSP) 这 类 离散 组 合 优化 问题 中 


迭代 次 数 


标准 测试 问题 LA16- 寻 优 曲 线 


二 总 总 四 名 
加 S 上 ba 
2 © S 


最 大 完工 时 间 (min) 


吉 
引 


1050 


1000 


950 
迁 代 次 数 


E 了 WOA 算法 解决 JSP 问题 的 可 行 性 ， 


部 分 标准 测试 问题 的 寻 优 


院 


并 


发 现 了 其 在 解决 JSP 问题 时 的 缺陷 。 为 了 克服 WOA 算法 解决 
JSP 问题 时 的 缺陷 


， 本 文 提出 了 


种 量子 鲸鱼 优化 算法 


迭代 次 数 


标准 测试 问题 LA31- 寻 优 由 线 


(QWOA), 并 对 其 进行 了 计算 复杂 度 分 析 、 收 合 怕 


仿真 分 析 。 
通过 对 
问题 时 ， 


迭代 次 数 


线 对 比 


点 


No 


在 未 来 


证 明 及 实验 


比分 析 发 现 , 本 文 所 提出 的 QWOA 算法 在 解决 JSP 
tk 有 和 迭代 次 数 少 、 收 敛 精 


全 局 搜索 能 力 强 等 优 


的 工作 中 ， 本 文 将 利用 WOA 算法 对 多 目标 组 合 优 


化 问题 进行 


究 ， 并 进一步 提升 WOA 算法 的 应 用 性 能 。 
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